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The ladder-rainbow truncation of the set of Dyson-Schwinger equations is used to study a variety of 
electroweak and strong processes involving light mesons. The parameters in the effective interaction are 
constrained by the chiral condensate and R; the current quark masses are fitted to TOtt and rriK. The 
obtained electromagnetic form factors are in good agreement with the data. Also the weak decay 
and the radiative and strong decays of the vector mesons agree reasonably well with the data. Finally, 
we indicate how processes such as tt-tt scattering can be described within this framework as well. 

1 Introduction 

Our goal is to describe the hadrons and their interactions in terms of their constituents, quarks 
and gluons, using the underlying theory, QCD. The set of Dyson-Schwinger equations [DSEs] 
form a useful tool for this purpose |j|]. In rainbow-ladder truncation, they have been successfully 
applied to calculate the masses and decay constants of light pseudoscalar and vector mesons [^,|3|. 
The dressed-quark propagator, as obtained from its DSE, together with the meson Bethe-Salpeter 
amplitude [BSA], form the necessary elements for calculations of strong interactions in impulse 
approximation, such as the p —> tttt decay. For electroweak processes, such as the electromagnetic 
form factors, radiative decays, and semileptonic decays, one also needs the qqy and qqW vertices. 

1.1 Dyson-Schwinger equations 

The DSE for the renormalized quark propagator in Euclidean space is 

S(p)~^ = iZ2^-\-Z4m(g)Zi J-^^g^D^,,{p-q)^y^S{q)ri{q,p) , ( 1 ) 

where is the dressed-gluon propagator and Pl,{q',p) the dressed-quark-gluon vertex. The 

most general solution of Eq. O has the form S{p) ^ = ifA{p^) + B{p^) and is renormalized at 
spacelike according to A{p^) = 1 and B{p?‘) = m{p) with m{p) the current quark mass. 

Mesons are described by solutions of the homogeneous BSE 

/ j4 

j^^K{p,q;Q) S{q+)TH{q+,q--,Q) 5’(g_), (2) 

at discrete values of = —najj, where mn is the meson mass. In this equation, p^. = p + pQ 
and p- = p — (1 — p)Q are the outgoing and incoming quark momenta respectively, and similarly 
for q±. The kernel K is the renormalized, amputated qq scattering kernel that is irreducible with 
respect to a pair of qq lines. Together with the canonical normalization condition for qq bound 
states, Eq. ( 0 ) completely determines the bound state BSA Tjy. Different types of mesons, such as 
(pseudo-)scalar, (axial-)vector, and tensor mesons, are characterized by different Dirac structures. 

The dressed qqy and qqW vertices satisfy an inhomogeneous BSE: e.g. the quark-photon vertex 
r^(p^,p_; Q), with Q the photon momentum and p± the quark momenta, satisfies Q 

r^(A+,P-;Q) =-^2 7 /^ + J-^i:^K{p,q-,Q) S{q+)Tf,{q+,q_-,Q) S{q-). (3) 

Solutions of the homogeneous version of Eq. (^) define vector meson bound states at timelike photon 
momenta = —rriy. It follows that r^(p+,p_) has poles at those locations 
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Table 1: Overview of results for the light pseudoscalar and vector meson masses and leptonic decay 
constants, all in GeV. Experimental data are from Ref. |^. _ 




fn 

rriK 

Ik 

rrip 

fp 

rriK* 

fx* 


f4> 

calc. 

0.138 

0.131 

0.497 

0.155 

0.742 

0.207 

0.936 

0.241 

1.072 

0.259 

expt. 

0.138 

0.131 

0.496 

0.160 

0.770 

0.216 

0.892 

0.225 

1.020 

0.236 


1.2 Model truncation 

To solve the BSE, we use a ladder truncation, 

K{p,q;P) -Q{{P - q?) , (4) 

in conjunction with the rainbow truncation for the quark DSE: T]^{q,p) —> 7i/A*/2 together with 
Zig^D^,{k) ^ in Eq. (§. This truncation preserves, independent of the details of the 

effective interaction both the vector Ward~Takahashi identity [WTI] for the qq'y vertex and 

the axial-vector WTI. The latter ensures the existence of massless pseudoscalar mesons connected 
with dynamical chiral symmetry breaking |^. In combination with impulse approximation, the 
former ensures electromagnetic current conservation |^. 

Eor the effective quark-antiquark interaction, we employ the Ansatz given in Ref. [^. The 
ultraviolet behavior of this effective interaction is chosen to be that of the QCD running coupling 
a(/c^); the ladder-rainbow truncation then generates the correct perturbative QCD structure of 
the DSE-BSE system of equations. In the infrared region, the interaction is sufficiently strong to 
produce a realistic value for the chiral condensate of about (240 GeV)^. With this model, we can 
solve the BSE for pseudoscalar and vector mesons, and calculate the meson masses and leptonic 
decay constants. The model parameters, along with the quark masses, are fitted to give a good 
description of the chiral condensate, and The results of our model calculations are 

shown in Table and are in reasonable agreement with the data. 

2 Meson interactions 

In impulse approximation, processes such as electromagnetic scattering, the weak Ki^ decay, radia¬ 
tive and strong decays of vector mesons, can all be described by the same generic loop integral 

r'^(P,Q,R) =iV, J^TT[Syq)T'^\q,q'-,P)S\q')T’’%q\q''-,Q)S<^{q'')T'^~^{q'^,q-,K)] , (5) 

where q —q' = P, q' — q" = Q, q" — q = K , and momentum conservation dictates P-|- Q -|- A' = 0. In 
Eq. d^, 5* is the dressed quark propagator with flavor index i, and P^{k, k'] P) stands for a generic 
vertex function with incoming quark flavor j and momentum k', and outgoing quark flavor i and 
momentum k. Depending on the specific process under consideration, this vertex function could 
be a meson BSA, a qqj vertex, or, in case of weak processes, a qqW vertex. In the calculations 
discussed below, the propagators, the meson BSAs, and the qq'y and qqW vertices are all obtained 
as solutions of their respective DSE in rainbow-ladder truncation, without adjusting any of the 
model parameters. 

2.1 Electromagnetic form factors 

Meson electromagnetic form factors in impulse approximation are described by two diagrams, with 
the photon coupled to the quark and to the antiquark respectively. Each diagram corresponds 
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Figure 1: On the left, our result for Q‘^F.,^{Q‘^), and right, our curve for the Ki^ form factor /+((5^). 


to an integral like Eq. @ with two meson BSAs and one gg 7 -vertex. With Q being the photon 
momentum, and the incoming and outgoing pseudoscalar mesons having momentum P =p Q/2, we 
can define a form factor for each of these diagrams |^] 

= lf^{P-Q/2,Q,-{P + Q/2)). (6) 


We work in the isospin symmetry limit, and thus F-n-iQ'^) = Fuuu{Q‘^)- The K~^ and form factors 
are given by = IFusu + IFuss and F^o = -^Fdsd + ^Fdss respectively. 

Our result for Q^F.,^ is shown in Fig. [l|, together with experimental data from Refs. the 

corresponding charge radius, together with the neutral and charged kaon charge radii, are given 
in Table The obtained charge radii agree quite well with the experimental data ||7, 10, as 
do our form factors. Up to about = 2GeV^, our result for F-j^ can be described very well by a 
monopole with mass scale given by our calculated nip, m = 742 MeV. Above this value, our curve 
starts to deviate more and more from this naive VMD monopole. Our results for Fk are given in 
Ref. p and can be fitted quite well up to about = 2 ~ 3 GeV^ by a monopole with mass scale 
slightly larger than the p mass. Asymptotically, these form factors behave like Q^F{Q'^) —> c up to 
logarithmic corrections |12|. However, numerical limitations prevent us from accurately determining 
the constants c. 


2.2 Weak interactions 

The matrix element {tt~{P + Q/2)\s'ypu\K^{P — Q/2)) describing the semileptonic decay of neutral 
kaons via a VF-boson with momentum Q can be characterized by two form factors 

/rf-(P_g/ 2 ,Q,-(P + Q/ 2 )) = 2PpUiQ^) + Qpf_iQ^) . (7) 


Table 2: Galculated charge radii in fm^, with expt. data and Ki 3 observables. The double 

entries for the expt. Ki^ data correspond to the neutral and charged Ki^ decays respectively. 



' TT 

^K+ 

rj^o 

A+ 

Aq 


r(iFe3) 

kiKps) 

calc. 

0.45 

0.38 

-0.086 

0.027 

0.018 

0.11 

7.38-10*^5-^ 

4.90-10*^5-^ 

expt. 

0.44 

0.34 

-0.054 

.0276, .0288 

.006, .025 

0.31, 0.11 

7.50, 3.89 

5.26, 2.57 
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The form factors f± for the decay are essentially the same as those for in the isospin limit, 
the only difference between the matrix elements for the and the K~^ decay is a factor of \/ 2 , 
the TT^ being {uu — dd)jV2, which results in a factor of 2 difference in the partial decay width. 

In the right panel of Fig. we show our result for /+(t) [^], together with the experimental 
data for this form factor |jl^. Experiments are often characterized in terms of the transverse, /+, 
and the scalar form factor /o, rather than /_, which is defined by 


/o(Q') = 


Q^7^^“(P-g/2,Q,-(P + g/2)) 


m 


K 


- mt 


= f+iQ") - 




m\ — ni^ 


f-{Q" 


( 8 ) 


The dimensionless slope parameter A for these form factors is defined as A = —m^/'(0)//(0); ^ 
is defined as /_(0)//+(0). The partial decay width can be obtained by integrating the decay rate, 
which depends on the lepton masses. Both the shape and the magnitude we obtain for these form 
factors agree well with experiments, as can be seen from Table 


2.3 Radiative decay of vector mesons 

We can describe the radiative decay of the vector mesons using the same loop integral, Eq. (I) , this 
time with one vector meson BSA, one pseudoscalar BSA, and one ( 7 ^ 7 -vertex []^| . The on-shell 
value gives us the coupling constant, which can be used to calculate the partial decay width. For 
virtual photons, we can define a form factor Fyp'yiQ'^), normalized to 1 at = 0 , which can be 
used in estimating meson-exchange contributions to hadronic processes IMl- 


In the isospin limit, the and vertices are identical, and are given by 


^ TUUU 

3 


{P,Q,-{P + Q)) = ^e^^o^pPaQpFpn^iQ^), 

rUn 


(9) 


where P is the p momentum. The a; vr y vertex is a factor of 3 larger, due to the difference in isospin 
factors. For the K* —> ATy decay, we have to add two terms: one with the photon coupled to the 
s-quark and one with the photon coupled to the u- or d-quark, corresponding to the charged or 
neutral K* decay respectively. 

As Eq. @ shows, it is gvP'y/mv that is the natural outcome of our calculations; therefore, it 
is this combination that we report in Table together with the corresponding decay widths |18|. 
The agreement between theory and experiment for gvp-tirnv is within about 10%, except for 
the discrepancy in the charged K* — > ATy decay for which we have no explanation. Likewise the 
large difference between the neutral and charged p decay width is beyond the reach of the isospin 
symmetric impulse approximation. Note that part of the difference between the experimental and 
calculated decay width comes from the phase space factor because our calculated vector meson 
masses deviate up to 5% from the physical masses. 


2.4 Strong decays of vector mesons 

If we continue the calculation of the electromagnetic form factors into the timelike region, we find a 
pole at the mass of the vector meson bound states. Using the behavior of the electromagnetic form 


Table 3: Vector meson radiative decays: coupling g/m in GeV ^ and partial decay width in keV. 



g/m 

rp=*=7r='='7 


g/m 

r 

UJTT'y 

g/m 

f'K*±K±-y 

g/m 

f'K*0K°'y 

calc. 

0.69 

53 

(53) 

2.07 

479 

0.99 

90 

1.19 

130 

expt. 

0.74 

68 

( 102 ) 

2.31 

717 

0.83 

50.3 

1.28 

116 
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Table 4: Overview of our results for vector meson strong decays (left): dimensionless coupling 
constants and partial decay width in MeV, and right, tt-tt scattering lengths, compared to leading 
order chiral perturbation theory, Oq = 7m‘^/{8irf^) and Oq = —m^/(47r/^) (Weinberg’s limit). 



dplTTZ 

r 

pTTTT 

gpKK 


gK*Kn 


calc. 

5.4 

115 

4.3 

6.7 

4.0 

31 

expt. 

6.02 

151 

4.64 

2.2 

4.60 

50 



Og Ug 

calc. 

Ref. [pq] 

0.170 0.045 

0.156 0.045 


factors Fuuu and around this pole, we can extract the coupling constants gpT^+-n- and g^px+K- 
respectively ||l^ , which govern the strong decays p tt-tt and (j) KK. The results from this 
analysis are given in Table ^ and are reasonably close to the experimental data. Similarly, the two 
form factors /+ and /o describing the weak Ki^ decay exhibit poles at = —rn?xi, and 
respectively, due to vector and scalar us bound states |j^. From the behavior close to the pole we 
can extract the coupling constant for the strong decay K* —> FT-tt as well. A direct calculation of 
the strong vector meson decays, using on-shell meson BSAs but different numerical techniques ||19|| , 
agrees reasonably well with these results extracted from the electroweak form factors. Note that 
the factor of three difference between the experimental and our calculated decay width for the (p is 
due to the phase factor (1 — 4m|^/m^)^/^: with our calculated masses, this factor is 0.051, whereas 
with the actual physical masses this factor is 0.015. The dimensionless coupling constants agree 
within 10% to 15% with the experimental data. 


3 TT-TT scattering 


Although impulse approximation seems to work remarkably well for a variety of interactions in¬ 
volving three external particles, one has to go beyond impulse approximation in order to describe 
processes with four (or more) external particles. As an example, consider tt-tt scattering at thresh¬ 
old. The generic loop integral for tt-tt scattering in impulse approximation is 

41 = w|^Tr[,s(g)r^(g,gO%')r.(9',g")%")r.(9",g'")^(r)r.(r,g)]. (lo) 


At threshold in the chiral limit, TT^{k + F/2, k — P/2) —> i'^^B{k ‘^)/and thus 


A 


4W 


/ 


d^k B\k)/ft 

(2Tr)4 (A:2A2(A;) +F2(A:))2 


( 11 ) 


which is nonzero. On the other hand, chiral symmetry dictates that the threshold scattering am¬ 
plitudes vanish in the chiral limit like //2 . Clearly, impulse approximation is insufficient to 

describe tt-tt scattering. 

In order to properly describe tt-tt scattering in the rainbow-ladder truncation, all possible 
diagrams with one or more insertions of the ladder kernel K across two pion BSAs should be 
added to the impulse contribution [^, as indicated in Fig. If we include these sets of ladder 
diagrams, we can show numerically that the threshold tt-tt scattering amplitudes indeed vanish like 
using the same model as in the previous section. The corresponding scattering lengths are 
given in Table We expect that in particular a)] will receive significant corrections from pion loop 
effects, which we have not included in our calculation: in chiral perturbation theory, higher order 
corrections (i.e. pion loops) change the leading order result to Ug = 0.220 and Og = 0.0444 

So far, we have only considered tt-tt scattering, but also in other hadronic 4-particle processes 
one should consider the contributions from these infinite sums of ladder terms, in addition to 


5 















Figure 2: Diagrams needed to correctly describe vr-vr scattering in rainbow-ladder truncation |21|. 


the impulse term. In general, we expect the role of these summed ladder contributions to be less 
important than in vr-vr scattering, except for processes that receive significant contributions from 
resonances. By adding these ladder diagrams one can unambiguously incorporate qq bound state 
effects, and we expect that this approach can provide a fundamental underpinning to many processes 
described by effective meson lagrangians. 
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